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ABSTRACT 

We report a new analysis of ASCA data on the iron K line complex in NGC 1068. 
The line complex basically consists of three components, as previously reported. A 
weak red wing of the 6.4 keV fluoresence iron K line is found. A plausible explanation 
is Compton scattering in optically thick, cold matter which can be identified with 
an obscuring torus or cold gas in the host galaxy. We also show that this 'Compton 
shoulder' should be observable with ASCA using a simulated reflection spectrum. In 
order to explain the two higher energy lines as well as the cold 6.4 keV line, we fit the 
ASCA data with a composite model of cold and warm reflection. This shows that cold 
reflection dominates the observed X-ray emission above 4 keV; the estimated scattering 
fraction from the warm medium is found to be ~ 0.08 per cent, an order of magnitude 
below the standard value of 1 per cent obtained from previous observations in other 
wavebands. The two higher energy lines have large equivalent width (~ 3 keV) with 
respect to the warm-scattered continuum, suggesting that efficient resonant scattering 
operates. The line energies are systematically lower than those expected from resonant 
lines for FeXXV and FeXXVI by ~ 100 eV. The redshifts may be due to either the 
ionized gas of the warm mirror receding at a radial velocity of 4000-5000 km s _1 , or 
effects of Compton scattering in a complicated geometry. 
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1 INTRODUCTION 

Since the discovery of the polarized broad line region (An- 
tonucci & Miller 1985), the standard model for the Seyfert 2 
galaxy NGC 1068 is of a Seyfert 1 nucleus hidden behind a 
thick obscuring torus. There is much observational support 
for this model (e.g., Antonucci 1993). Contrary to the pow- 
erful nuclear activity observed at some other wavelengths, 
the weakness of the observed X-ray luminosity and failure 
to detect any direct component, even up to 20 keV, strongly 
suggests that the direct X-ray radiation from the central 
source is blocked entirely from our direct view (Koyama et 
al 1989; Mulchaey et al 1995). 

The lack of low energy absorption in the X-ray spec- 
trum was interpreted as evidence that the observed X-rays 
are the scattered radiation of a hidden Seyfert 1 nucleus 
(e.g., Elvis & Lawrence 1988). However, extended soft X- 
ray emission (~ 13 kpc) is resolved with the ROSAT HRI 
(Wilson et al 1992), and an ASCA spectrum shows ther- 
mal emission, probably coming from the spatially extended 
circumnuclear starburst, dominates in the soft X-ray band 
(Ueno et al 1994). The ASCA spectrum then no longer looks 



like a Seyfert 1, indicating very steep thermal emission below 
3 keV (which now agrees with the Einstein Observatory IPC 
measurement by Monier & Halpern 1984) and only above 3 
keV, where a flat continuum (r ~ 1.2) with strong iron K 
lines (Ueno et al 1994) is seen, is emission presumed to orig- 
inate from the active nucleus. 

Based on the unification model, a large equivalent width 
(EW) of the iron K line was predicted by Krolik & Kall- 
man (1987). This was later verified by a Ginga observation 
(EW ~ 1.3 keV, Koyama et al 1989; Smith et al 1993). The 
line centroid energy of the iron K line was 6.55 ±0.1 keV 
which suggested a blend of cold and ionized lines, but it 
could not be resolved at the spectral resolution of the Ginga 
LAC (Turner et al 1989). A subsequent BBXRT observation 
resolved the iron K line feature into three components for 
the first time (Marshall et al 1993). They were identified 
with cold Fe (less ionized than FeXVII at 6.4 keV), FeXXV 
and FeXXVI. 

Several attempts to explain the iron K line feature have 
been made (Band et al 1991; Marshall et al 1994), assuming 
that all X-ray emission is scattered light from the active nu- 
cleus and the scattering medium is optically thin in order to 
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Table 1. Four gaussian fit to the iron K line complex. The line 
width of each component is assumed to be u = 10 eV, since 
no significant broadncning (<r < 80 eV) is found. The EWs are 
with respect to the single power-law continuum with Galactic 
absorption (see text). Line-intensity ratios to the 6.4 keV line 
(I/I6.4) are also shown. 



E 
keV 

6.21 ±0.16 
6.40 ±0.03 
6.61 ±0.03 
6.86 ± 0.05 



I 

10 _5 phs _1 cm" 



EW 
keV 



I/I 6 .4 



0.51 
5.19 
3.20 
2.65 



+0.70 

0.50 
+0.79 
1.27 
+ 1.00 
0.80 
+0.89 
0.87 



0.11 
1.10 

0.69 
0.58 



0.15 
0.10 
+0.17 
0.27 
0.22 
0.17 
+0.19 
0.19 



0.10 



+0.13 



0.09 
1.00 



0.62 
0.51 



+0.23 
0.18 

+0.20 
0.20 



account for the lack of low energy absorption. The observed 
continuum above 3 keV is unusually flat compared with that 
in Seyfert 1 galaxies observed with ASCA (F ~ 1.9; Nandra 
et al 1996; Reynolds 1996). Together with the presence of a 
strong 6.4 keV line, this instead suggests that reflection from 
cold, thick matter is an important component, as observed 
in other Seyfert 2 galaxies such as NGC 6552 (Fukazawa et 
al 1994) and the Circinus galaxy (Matt et al 1996a) . To ac- 
count for the three line components, Matt, Brandt & Fabian 
(1996) proposed a composite scenario in which two distinct 
reflections from an optically thick torus and an optically 
thin, warm scattering mirror together compose the observed 
X-ray spectrum of NGC 1068. 

We report here the discovery of a Compton shoulder on 
the 6.4 keV fluorescence iron K line and a spectral analysis 
of the ASCA PV data based on the model by Matt et al 
(1996b). We assume the distance of NGC 1068 to be 22 Mpc 
throughout this paper. 



2 THE ASCA DATA 

NGC 1068 was observed on 1993 July 24 with ASCA 
(Tanaka, Inoue & Holt 1994) during the Performance Veri- 
fication (PV) phase. Initial results have been published by 
Ueno et al (1994). We use only the Solid state Imaging Spec- 
trometer (SIS; SO and SI) data, because of their superior 
spectral resolution for resolving the iron K line complex. 
The source data of 38 ks are extracted from a circular re- 
gion of a radius 3.2 arcmin, and the background data are 
taken from a source-free field in the same detector. Spec- 
tral fits are performed with the SO and SI data jointly while 
the figures presented are made by summing both detectors 
for display purposes. Errors on each spectral parameter are 
given at the 90 per cent confidence level for one interesting 
parameter. 



2.1 Compton shoulder 

We show the observed line profile of the iron K complex 
in Fig. 1. The underlying continuum is fitted with a sin- 
gle power-law model described below. Since X-ray emission 
below 3 keV is dominated by starburst emission (Ueno et 
al 1994) and there is evidence for emission lines due to Ar 
and Ca (see below) between 3 and 4 keV, we use the 4- 
10 keV data. Fitting the 4-10 keV data excluding the iron 
K line band 6-7 keV with a single power-law plus Galactic 
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Figure 1. (a) Total line profile of the iron K line complex in 
NGC 1068 obtained from the ASCA SIS. The SO and SI arc 
summed to display this figure. Note that the instrumental effi- 
ciency declines towards higher energy; e.g., the efficiency at 6.9 
keV is 15 per cent lower than at 6.4 keV. (6) The two higher 
energy lines are subtracted from the total line profile. A clear 
line asymmetry is seen. The red wing can be identified with a 
Compton shoulder of the 6.4 keV line. 



absorption Ah— 4.5 x 10 20 cm" 2 (Stark et al 1992), leaves 
the total line profile shown in Fig la. The photon-index is 
r = 0.35t o ^, flatter than that given in Ueno et al (1994) 
because of the exclusion of the 3-4 keV band where some 
emission lines are present, causing the slope to steepen. 

As previously suggested by Marshall et al (1993) and 
Ueno et al (1994), three emission-line peaks are seen. These 
three lines do not show evidence for significant broadening 
(a < 80 eV). However, there is a weak wing on the red side 
of the 6.4 keV line. To clarify this, we show only the 6.4 keV 
line profile in Fig. lb, which was produced by fitting the two 
higher energy lines with a narrow (er = 10 eV) gaussian and 
subtracting them from the total line profile. The subtrac- 
tion of the two lines is fairly good as there are no significant 
residuals seen above 6.5 keV. The line profile is clearly more 
complicated than a single gaussian; the core emission is con- 
sistent with the SIS energy resolution, whereas a wing can 
be seen on the red side in contrast. Fitting the red wing 
by a gaussian gives a line energy of 6.21 ± 0.16 keV and an 
intensity one tenth of the 6.4 keV line with ~ 90 per cent 
confidence level. The four gaussians plus a power-law model 
now gives a good fit to the 4-10 keV data with \ 2 = 54.98 for 
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Figure 2. (a) Line profile of an iron K line from a Monte Carlo 
simulation of reflection in thick, cold matter (see text). To show 
the shape of the weak Compton shoulder, the intensity here has a 
logarithmic scale, (b) The line profile (solid line) folded through 
the SIS response from (a). A narrow (o = 10 eV) gaussian with 
an intensity normalized to the peak intensity of the line core is 
also shown by the dotted line. Note that the spectral resolution of 
the SIS is FWHM 130 eV at 6.4 keV, appropriate for our data 
which were taken during the PV phase, (c) Difference between the 
line with the Compton shoulder and a narrow gaussian in (b). It 
is clear that the residual is observable with the SIS, consistent 
with the observed result. 



70 degrees of freedom. The line energies, equivalent widths 
and relative line intensities of each line component are sum- 
marized in Table 1. 

The possibility of a symmetric broad base to the 6.4 
keV line cannot be excluded, since the subtraction of the 
higher energy lines will have eliminated any blue wing. A 
double gaussian, with a narrow and a broad component, fit 
to the 6.4 keV line profile is slightly worse than the one with 
the red wing (by A\ 2 — 1-4), but acceptable. In this case, 
the line width of the broad component is a ~ 0.3 keV. 

Compton scattering is a likely mechanism to explain 
the red wing. A line extending only to lower energy indi- 



cates that the electrons responsible for scattering are cold 
(Pozdnyakov, Sobol & Sunyaev 1979). Theoretical calcula- 
tions of Compton reflection spectrum from optically thick, 
cold matter have already predicted such a weak red wing, 
the so called 'Compton shoulder' (George & Fabian 1991; 
Reynolds et al 1994), which is a characteristic feature due to 
Compton scattering of the 6.4 keV fluorescent line photons. 
It consists of a series of shoulders extending towards low en- 
ergies. Computed line shapes due to the Compton scattering 
are presented in Hatchett & Weaver (1977); Pozdnyakov et 
al (1979); George & Fabian (1991); and Sunyaev & Chu- 
razov (1996). The first shoulder is due to single scattering 
and is spread over two Compton wavelengths (~ 160 eV); 
multiple scattering makes a further weak low energy tail. 

If pure reflection from thick, cold matter is observed 
with ASCA, the Compton shoulder should be resolvable 
with the SIS. To demonstrate this we have used a Monte 
Carlo code kindly provided by C. Reynolds (see also George 
& Fabian 1991). The line profile derived from a Monte Carlo 
reflection spectrum with spectral resolution of 10 eV (Fig. 
2a) is compared, after folding through the SIS response, with 
that of a narrow (a — 10 eV) gaussian (Fig. 2b). The differ- 
ence bewtween the two profiles should be detectable (Fig. 
2c). The shape of the Compton shoulder depends slightly 
on geometry (e.g. Matt, Fabian & Reynolds 1996) and its 
height also depends on the abundance ratio of iron and the 
lighter elements, mainly oxygen, which absorb the 6.4 keV 
photons (George & Fabian 1991). In order that the shoulder 
is maximally broad, a large back-scattered component is re- 
quired. The quality of the present data is not good enough 
for detailed spectral modelling, however, and we merely as- 
sume that reflection occurs from a face-on, semi-infinite slab 
with an optical depth of unity. The strength of the simula- 
tion has been adjusted to the observed intensity level. Al- 
though the residual profile due to the Compton shoulder 
shows some asymmetry, fitting this with a gaussian gives 
the line centroid energy ~ 6.24 keV and the intensity ratio 
to the gaussian line-core ~ 12 per cent. These values are in 
good agreement with the gaussian modelling of the observed 
line profile (see Table 1). 

Note that the broad line from a relativistic accretion 
disk, such as that observed in Seyfert 1 galaxies (e.g., Tanaka 
et al 1995), is an unlikely explanation for the red wing, since 
we do not directly see the relativistic line emitting region 
(which emerges from within ~ 40 gravitational radii). A 
scattered relativistic line is also implausible. The scattering 
cone appears to have a half opening angle of ~ 40° (Evans 
et al 1991; Tsvetanov et al 1996) and the inclination angle of 
the accretion disk in NGC 1068 with respect to most of the 
scattering volume must therefore be near 40° . In this geom- 
etry the red wing of the diskline would be skewed down to 4 
keV and peaked at energies well below 6 keV (see inclination 
angle dependence of diskline profiles in Fabian et al 1989) 
giving rise to a much broader red wing than observed. 

2.2 Spectral fit with the composite model 

The presence of the Compton shoulder indicates that reflec- 
tion from optically thick cold matter is an important process 
in producing the observed hard X-ray emission in NGC 1068. 
However, pure cold reflection alone does not explain the two 
higher energy lines. 
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Figure 3. ASCA SIS spectrum of NGC 1068 fitted with the com- 
posite model. The best-fit model for the 3.2-10 keV band is shown 
as well as the cold reflection and warm scattered continuum plus 
ionized lines. The data below 3.2 keV down to 2.3 keV are also 
shown so that excess emission due to the starburst can be seen. 




(cold) [10 40 erg/s] 



Figure 4. Contour plot of the 2-10 keV luminosities of the cold 
and the warm reflection continua obtained from the composite 
model fit. The luminosities are derived assuming the distance of 
NGC 1068 to be 22 Mpc. Note that the cold reflection luminosity 
is integrated over 27r (see text). Contours arc drawn at 68 and 
90 per cent confidence levels for two interesting parameters. The 
best-fit values, 2L C = 1.36 X 10 41 crg s" 1 and L w = 5.47 X 
10 40 crg s _1 , arc shown by a plus symbol. 



Table 2. Equivalent widths of two calcium-K and three iron- 
K lines respect to each continuum conponent obtained from the 
composite model fit. EWq and EWw are equivalent widths re- 
spect to the cold reflection continuum and the power-law warm- 
reflection continuum, respectively. Note that the EW for the 6.4 
keV line is the sum of the line-core and Compton shoulder. 



Eline 

keV 

3.71 
3.94 
6.40 
6.61 
6.86 



EW C 
keV 



+0.11 



0.19 



i 91 +0.26 
^ -0.28 



EW W 
keV 



11+ 006 
u - ±1 -o.io 



J -^ y -0.87 
9 q 9 +0.98 



It is likely that the higher energy lines are produced in 
optically thin ionized (warm) matter which is located above 
the obscuring torus and electron-scatters the continuum of 
the central source, as originally proposed by Krolik & Kail- 
man (1987). If the ionized scattering medium is optically 
thin enough, photoelectric absorption has little effect on the 
electron-scattered continuum so the emergent spectrum re- 
tains almost all the original shape of the central source, in 
the ASCA band. 

Thus, as proposed by Matt et al (1996b), a combination 
of cold and warm reflections could be a plausible description 
for the observed hard X-ray emission from NGC 1068. We 
therefore fit the ASCA spectrum with the composite model 
described above. 

We assume that the incident spectrum of the central 
source is a power- law with V — 2.0, typical of Seyfert Is (e.g., 
Nandra et al 1996; Reynolds 1996). The pure cold reflection 
continuum and the warm-scattered continuum are approx- 
imated by the viewing-angle averaged Lightman & White 
(1988) type reflection model and a power-law of r = 2.0 
with no intrinsic absorption, respectively. The iron K line 
complex is modelled by four gaussians including the Comp- 
ton shoulder. Weak line features due to calcium between 
3.6-4 keV are also modelled by narrow gaussians. The 3.2- 
10 keV data are fitted with the composite model, and the 
best fit model gives x 2 = 70.14 for 83 degrees of freedom. 
The EW of each line are shown in Table 2. Note that the 
EW of the cold iron K line is the sum of the 6.4 keV line- 
core and the Compton shoulder measured with respect to 
the cold reflection continuum, and those of the two higher 
energy lines are with respect to the warm-scattered T = 2.0 
power-law continuum. 

Although only the 3.2-10 keV data are fitted, to show 
excess emission due to the starburst below the low energy 
bound, the 2.3-10 keV data are shown in Fig. 3. Two sig- 
nificant lines at 3.71 ± 0.04 keV and 3.94 ± 0.09 keV are 
identified with a neutral and helium- like calcium (CaXIX). 
The EW of the neutral Ca line with respect to the cold 
reflection continuum is EW = 190 ± 110 eV, which is sim- 
ilar to that in NGC 6552 (Fukazawa et al 1994) but larger 
than that expected from a Monte Carlo simulation (~ 50 
eV; Reynolds et al 1994). This is further evidence for cold 
reflection. Another line feature around 3.1 keV is due to 
Ka emission from ArXVII which is probably from both the 
thermal emission and the warm scattering matter. 

The 2-10 keV fluxes of the two components esti- 
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mated from the best-fit model are Fx (cold) = 2.5 x 
10~ 12 erg cm -2 s _1 for the cold reflection, and Fx(warm) = 
1.0 x 10~ 12 erg cm~ 2 s _1 for the warm scattered continuum. 
Confidence contours of the corresponding luminosities are 
given in Fig. 4, using a distance of 22 Mpc. The best-fit val- 
ues m the 2-10 keV band are 6.8 x 10 40 erg s" 1 for the cold 
reflection continuum and 5.5 x 10 40 erg s _1 for the warm re- 
flection continuum, respectively. It should be noted that the 
luminosity is integrated over 2tt for the cold reflection be- 
cause of the assumed geometry whilst over 4n for the warm 
reflection. However, for a direct comparison of observed in- 
tensity between the two components, the cold reflection lu- 
minosity multiplied by 2 is shown in Fig. 4. 

If significantly large absorption (say, Nn> 10 cm~ ) 
or a starburst component is included in the fit, the value of 
Lw would be reduced slightly. Fitting the data in a higher 
energy band (e.g. 4-10 keV) gives poorer constraints on Lc 
and Lw- Results are consistent with those in Fig. 4, although 
the warm reflection component is not necessary at the 90 per 
cent confidence level. 



3 DISCUSSION 

3.1 Intrinsic luminosity 

A rough estimate of the intrinsic X-ray luminosity can 
be obtained using the warm scattering luminosity Lw ~ 
5.5 x 10 40 erg s _1 derived from the composite model fit. For 
highly-ionized iron such as FeXXV and FeXXVI, resonant 
scattering plays an important role in making the EW large, 
if the scattering medium is optically thin (Band et al 1991; 
Matt et al 1996), but it is easily suppressed below the large 
observationally inferred values (Table 2) when the column 
density is, say, larger than 10 22 cm~ 2 (e.g., see Matt et al 
1996b). This indicates that the optical depth of the warm 
scattering mirror is small enough to allow efficient reso- 
nant scattering. Thus assuming a scattering depth (nZux) of 
0.002t .oo2 so that the column density of the scattering mir- 
ror is Nn~ 3x 10 21 T .oo2cm~ 2 , and the fraction of the sky oc- 
cupied by the scattering mirror, 0.25(fi/47r)o.25 , the intrin- 
sic luminosity is estimated to be L' = Lwt^ 1 (Q,/ Atv)^ 1 ~ 
1 x 10 44 r ( ^ 1 02 (f2/47r)Q 2 5 erg s _1 in the 2-10 keV band, where 
the density and depth of the scattering medium are n and I, 
respectively, and Lw is the measured luminosity of the warm 
reflection continuum. This estimate is entirely independent 
of data in other wavebands. 

There are further possible ways to estimate the intrin- 
sic luminosity. Using the correlation between X-ray and 
[OIII]A5007 luminosities (logL[ iii]/logL 2 -iokeV ~ —1.8) 
derived by Mulchaey et al (1995), the 2-10 keV luminos- 
ity would be V ~ 7 x 10 43 erg s" 1 . Note that the [OIII] flux 
was measured with a large aparture (10-15 arcsec, see Whit- 
tle 1992) containing the strong starburst region which could 
contribute some fraction of the flux. The estimate from the 
above, X-ray based, argument is consistent with this value. 
For a comparison, the far infrared luminosity measured with 
IRAS is 5 x 10 44 erg s _1 (e.g., David, Jones & Forman 1993), 
although this does contain a significant contribution from 
the circumnuclear starburst. 



3.2 Reflection from cold matter 

The detection of the Compton shoulder of the 6.4 keV iron 
fluorescent line is strong evidence that we are seeing reflec- 
tion from optically thick, cold matter. The presence of the 
neutral calcium line supports this. The flat 4-10 keV con- 
tinuum (r ~ 0.4) suggests that cold reflection is a major 
component of the observed X-ray emission from the active 
nucleus which dominates the observed radiation only in the 
energy band above 3 keV. In fact, the fit with the composite 
model of cold and warm reflections (Section 2.2) indicates 
that ~ 71 per cent of the total 2-10 keV flux comes from 
cold reflection. 

Presumably the cold reflection occurs at the inner sur- 
face of the optically-thick, obscuring torus which entirely 
blocks the direct radiation from the central source in our di- 
rection (Ghisellini, Haardt & Matt 1993; Matt et al 1996b). 
For such a geometry only part of the reflection surface is 
visible to us. If that fraction is fc, the observed luminosity 
is described as Lc — fcr/L' , where rj is the total albedo of a 
cold reflector subtending a solid angle of 2ir. Using values of 
L c = 6.8 x 10 40 erg s _1 (see Section 2.2 and Fig. 4), r? = 0.022 
(in the 2-10 keV band), and L 1 = 7 x 10 43 Li 3 . 8 erg s" 1 (see 
Section 3.1), fc is then found to be « 4.4 x 10 _2 L^_ 1 8 . 

The small value for fc suggests a high inclination of the 
obscuring torus (close to edge-on) if its inner surface is the 
reflecting source. Alternatively a torus of a small radial/axial 
ratio (Pier & Krolik 1993) is also possible, if it is not nearly 
edge-on. Recent water maser observations by Gallimore et al 
(1996) and Greenhill et al (1996) suggest a high inclination 
> 82 degrees of the maser disk which is on a pc scale with 
a toroidal geometry and is identified with the X-ray absorb- 
ing gas torus at large optical depth (i.e., Nn> 10 25 cm~ 2 , 
Koyama et al 1989). A rather broad [OIII]A5007 emission- 
line cone (projected opening angle plo j ~ 65-80°, Cecil, 
Bland & Tully 1990; Evans et al 1991) also favours an edge- 
on torus over one with a large scale-height. 

Another possible cold reflector is the cold gas in the 
bulge and disk of the host galaxy (Fabian 1977). Strong 6.4 
keV line emission is observed from molecular clouds in our 
Galactic centre, where a strong iron K line emitter is Sgr B2, 
~ 100 pc away from the Galactic centre. The X-ray spec- 
trum there also shows a complex of iron K lines (Koyama et 
al 1996). As the stellar disk of NGC 1068 is almost face-on, 
contrary to the torus, radiation from the central source es- 
caping from the torus can pass through some of the cold gas 
clouds in the host galaxy. The likely distances of cold clouds 
from the active nucleus and the disk-like geometry make 
the covering factor small, and consistent with the small fc. 
If this is the case, future high spatial resolution mission like 
AXAF would resolve a conical shape to the iron K line emis- 
sion. 

The cold reflection spectrum has been studied by many 
authors (Guilbert & Rees 1988; Lightman & White 1988; 
George & Fabian 1991; Matt, Perola & Piro 1991; Reynolds 
et al 1994). The cold reflection-dominated spectrum in 
NGC 1068 predicts that a spectral break due to Compton 
down-scattering would be observed around 30 keV. This is 
consistent with the CGRO/OSSE upper limit and may be 
observable with high quality data from RXTE and SAX. 
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3.3 Reflection from warm matter 

Cold reflection is a unique phenomenon at X-ray wave- 
lengths. Warm scattering is much more relevant to the ob- 
served optical/UV scattering. The 2-10 keV luminosity of 
the warm-scattered continuum is Lw ~ 5.5 x 10 40 erg s~\ 
which implies a scattering fraction in the warm matter to 
be /w = Lyf/L' ~ 7.8 x W~ 4 L'^g. Various other obser- 
vations of NGC 1068 suggest that the scattering fraction is 
approximately 1 per cent (summarized in Pier et al 1994). 
The X-ray scattering fraction appears to be an order of mag- 
nitude lower. Pier et al (1994) deduced the bolometric lumi- 
nosity of NGC 1068 to be - 7 x lO^/^erg s _1 , where the 
scattering fraction is 0.01 /o.oi- Since they assumed that the 
whole hard X-ray emission above 2 keV (as observed with 
BBXRT, Marshall et al 1993) is scattered by the warm mir- 
ror, contrary to our result, their estimate has been modified 
slightly. 

If the two higher energy lines are due to FeXXV and 
FeXXVI, the ionization parameter of the warm scattering 
matter, £ = L' /nR , is about several hundreds (Kallman 
& McCray 1982; Krolik & Kallman 1987; Band et al 1991), 
where R is a distance of the scattering medium from the 
central source. Assuming I = R, the column density can be 
written as nl = L'(£Z) _1 . Since /w = <7TnZ(f2/47r), 



If £ = 600^600 and fi/47r = 0.25(n/47r) .25 (e.g., Evans 
ct al 1991) are taken, then the size of the scattering re- 
gion is approximately I ~ 8£,qq Q L2 3 . 8 (£I/4it)o.25 pc. This 
is large enough to be visible to us without occultation by 
a pc-scale torus; direct X-ray imaging of this is impossi- 
ble at present as the angular size is the order of 0.1 arc- 
sec (» 11 pc at D = 22 Mpc). The size of the optical/UV 
scattering nebula observed by spectropolarimetry with the 
Hubble Space Telescope (HST) is several tens pc or up 
to 200 pc (Antonucci et al 1994; Capetti et al 1995), the 
warm X-ray mirror is therefore smaller. The column density 
nl — 5 x 10 21 I/43 1 8 (n/47r)o.25cm -2 , which is small enough to 
allow efficient production of the resonant scattering lines of 
FeXXV and FeXXVI as mentioned in Section 3.1. To main- 
tain the high value of £ in order to emit the highly ionized 
iron lines over a large scale, a low mean electron density such 
as ~ 200 cm~ 3 is required. 

3.4 Complexities 

We have discussed X-ray scattering from a warm medium in 
the context of the standard unification model as proposed 
by Krolik & Kallman (1987) and Matt et al (1996b), i.e., 
the scattering medium is located outside the obscuring torus 
and is directly visible to us. However, there are some prob- 
lems: (1) a very small scattering fraction /w ~ 0.8 x 10~ 3 ; 
(2) disagreement of line energies discussed below; and (3) 
the location of the region where the higher energy lines are 
formed. 

We note that the line energies for the two higher energy 
lines (6.61 keV and 6.86 keV) disagree with the energies of 
the resonant lines of FeXXV, 6.70 keV, and FeXXVI, 6.97 
keV, (see Table 1) whilst that of the cold iron K line agrees 
with 6.40 keV perfectly. The shift of the line energies are 



-90 ± 30 eV for FeXXV and -110 ± 50 eV for FeXXVI. We 
have checked that there are no significantly strong iron lines 
at the observed energies from any other ionization states 
(e.g., Bciersdorfer et al 1993). Interestingly, the shift of the 
two line energies are similar to each other (~ 1.5 ± 0.6 per 
cent lower than rest energies). If the redshift is due to the 
Doppler effect from receding warm scattering matter, then 
the radial velocity is ~ 4000-5000 km s _1 . The optical po- 
larized broad H/3 line is also redshifted but only by 600 
km s _1 (Antonucci &: Miller 1985) from the narrow line 
component, which is explained as scattering in outflowing 
wind (Krolik & Begelman 1986; Balsara & Krolik 1993). As 
the discrepancy between the velocities is too large, the same 
outflowing wind is perhaps not an explanation. It should be 
noted however that the X-ray lines are produced in the mir- 
ror itself whilst the polarized broad line is scattered by the 
mirror. 

The location of the gas emitting the two higher energy 
lines is an issue. It is supposed to lie outside the torus. Sys- 
tematic studies of Seyfert 1 galaxies with ASCA show that 
a large fraction of them show evidence for partially ionized 
gas in the line of sight, implied by absorption edge feature 
mainly due to oxygen, OVII and OVIII (Reynolds 1996). A 
recent study of the warm absorber in MCG-6-30-15 (Otani 
et al 1996) suggests that the OVII and OVIII absorbers oc- 
cupy distinct regions; the variability study locates the OVII 
absorber > 1 pc far from the central source, and the OVIII 
absorber in the broad-line region. If we observed this ob- 
ject from a different viewing angle, the OVII warm absorber 
would be seen as a scattering mirror. However, the ionization 
state implied from OVII (£ < a few tens erg cm s _1 ) is too 
low for the emission of FeXXV and FeXXVI. If ionized gas 
in the nucleus of NGC 1068 has a similar ionization struc- 
ture to that of MCG-6-30-15, the gas emitting the highly- 
ionized iron line should be located well within the torus like 
the OVIII absorber so that it would not be visible to us 
directly. If this is the case, the radiation emitted from the 
inner hot gas, including resonantly scattered FeXXV and 
FeXXVI lines, is also scattered into our line of sight by the 
cold reflecting matter. The line emission would be Compton 
scattered and then end up having line energies shifted lower 
as observed. Scattered continuum radiation would have the 
same shape as the cold reflection continuum. 

We do however know that the optical/UV mirror (see 
e.g., Miller, Goodrich & Mathews 1991) is located outside 
the torus and is at least directly visible to us, so we should 
see some X-ray continuum scattered from it. (It is unlikely 
that this mirror is sufficiently ionized to produce the high 
ionization X-ray iron lines.) This may reconcile the small 
value of fw- Importantly, the EW with respect to the warm 
reflection continuum, discussed before (e.g., Table 2, Sec- 
tions 2.2 and 3.1), is now meaningless since the spatial or- 
gin of the lines and continuum is not the same. Uncertainties 
in the geometry of the inner nucleus and the mirrors make 
it difficult to estimate the expected luminosities from each 
component. 

A further possibility, which can increase the predicted 
EW of the lines from highly ionized iron, is that there is 
velocity shear in the warm mirror. As discussed in Section 
3.1, the EW decreases as the column density of the mirror 
increases as the line is absorbed. This assumes that the scat- 
tering medium is all at the same velocity. If, however there is 



© 0000 RAS, MNRAS 000, 000-000 



Iron K line complex in NGC 1068 7 



shear such that the velocity doppler-shifts the line by more 
than its spectral width (natural and thermal) over a typical 
mean free path, corresponding to a column of ~ 10 21 cm~ 2 , 
then much of the line is not resonantly absorbed in other 
parts of the mirror. A larger EW then arises from greater 
optical depths than assumed in Section 3.1 (see Fig. 9 in 
Matt et al 1996 for an indication of the effect). For a mirror 
gas temperature of about 10 6 K, the velocity shear needs to 
exceed about 50km s _1 (10 21 cm 2 ) -1 for this to be signifi- 
cant. The upper limit on the width of the line (< 80 eV; 
Table 1) limits the extent of this effect. In this case the in- 
trinsic luminosity can be less than we estimated in Section 
3.1. 

Higher resolution spectroscopy resolving the individual 
line shapes will test and discriminate between these hypoth- 
esis. Given the complex geometry and the possible complex- 
ities in the formation of the iron and other X-ray lines, we 
expect it to reveal that NGC 1068 has a rich and interesting 
spectrum. 
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